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ABSTRACT 

Measurements of the microlensing optical depth t towards the Galactic bulge 
appear to depend on the method used to obtain them. Those values based on the 
lensing of red clump giants (RCGs) appear to be significantly lower than those based 
on the lensing of all stars along the line of sight. This discrepancy is still not under- 
stood. Through Monte Carlo simulations, it is found that the discrepancy cannot be 
explained by a dependance on the flux limits of the two methods. The optical depth 
is expected to be generally constant as a function of source apparent magnitude for 
Io > 13.0, except in the range 13.5 < Iq < 15.5. Here many RCGs are detected, 
causing a significant oscillation in r. The amplitude of this oscillation is a function 
of the inclination angle of the Galactic bar, (9b ar , which may thus be constrained. A 
further constraint comes from a similar dependance of r with dbuY- combining the 
predicted trends with the measured values provides ler upper limits, which exclude 
the large bar angles recently reported by the GLIMPSE and EROS surveys. The 
latest survey data from EROS-2 appear to show the predicted t oscillation, though 
currently at a low significance. However, a further sign comes from EROS-2 event 
counts, which show a clear skew towards fainter magnitudes. 
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1 INTRODUCTION 

Thousands of microlensing even ts in our Galaxy have so far 



CHO (e.g. 



and EROS (e.g. 



Thomas et alj 


2005 


), MOA (e.g. 


Sumi et al. 


2003) 



Afonso et alj2003T) . with many more detected 



every yeaj^. One of the most important measurements that can 
be made from these observations is of the optical depth, r - the 
probability of seeing a microlensing event. However, the mea- 
sured value appears to depend strongly on the method used to 
obtain it. 



1 Real-time alerts are available online, e.g. the OGLE Early Warning 
Email: A.A.Wood@postgrad.manchester.ac.uk System: www.astrouw.edu.pl/~ogle/ogle3/ews/ews.html 
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Popowski et al. 



2003) reported tmacho = 2.17±|£ 



10" 6 at (l,b) = (1. 50°, -2.68° ) , and more recently from 
the OGLE-II survey, 



Sumi et al 



1 2006) found togle = 
2.55ig;^ x 10" 6 at (l,b) = (1.16°, -2.75°). The MA- 
CHO value was based on the lensing of 42 red clump giants 
(RCGs), and used standard photometric fitting. The OGLE 
analysis instead used difference image analysis (DIA), but 
was similarly based on the lensing of 32 red giants, red su- 
per giants and RCGs, and obtained an optical depth consis- 
tent with the previous MACHO result. However, both these 
values are significantly lower than two other recent mea- 
surements, which were based on t he len sing of all stars. 
Using 28 DIA events, 



Sumi ct al 



1 2003T) found tmoa 



3.361*;^ x 10~ 6 [0.77/(1 - / disc )], where / disc is the con- 
tribution from disc sources - the coordinates of this value 



are given in 



Sumi et al 



Alcock et al 



1 20061) : (l,b) = (3.0°, -3.8°). 
1 20001) had previously found r = 3.23±g j§ x 
10~ 6 [0.75/(1 - /disc)] at (l,b) = (2.68°, -3.35°), from 
99 MACHO DIA events. The latest measurement comes from 
the EROS-2 survey of b ulge RCGs, which yielded 120 events: 

= (1.62 ± 

3 deg)] x 10~ 6 , where a = (0.43 ± 



Hamadache et al 



20061) give the trend teros 

0.23) exp[-a(]&] 
0.16 deg- 1 ), in the latitude range 1.4° < |fe| < 7.0°. This 
agrees well with previous EROS values, and with the recent 
MACHO and OGLE-II measurements. 

The question naturally arises as to why the RCG-based 
optical depths appear to be lower than those from all stars. One 
possibility is a dependence on the flux limits of the two meth- 
ods. RCGs are bright; the latter method will include much 

fainter stars, and so probe sources at greater distances, which 

I | I 
will have a higher optical depth (Stanek 1995). 

This potential explanation of the discrepancy is investi- 
gated using Monte Carlo simulations of Galactic microlens- 
ing, and the optical depth as a function of source apparent 
magnitude is then predicted. The model Galaxy is barred, and 
in light of observations by the Spitzer Space Telescope (SST), 
that support a bar inclination angle much larger than suggested 
by previous studies (see ^24}, the effect on the expected r of 
changing the bar angle is determined. Combining these results 
with the observed optical depths, upper limits are placed on 
the bar angle. Sj2]describes the model. The results and discus- 
sion are presented in Sj3] the model results are given in i )3.1| 



and comparisons are made with the latest EROS-2 data in $3.2\ 
A summary and conclusions follow in <J4] 



2 THE MODEL 

2.1 Bulge and disc mass models 

The mass models and param eters of the Galactic bulge (bar) 



and disc are as described in 
based on those of 



normalised the G2 bulge model of 



Han & Gould ( 2003), who emp irically 



Wood & Mao (2005). They are 



Dwek et al 



1995, table 



1) with Hubble Space Te lescope star counts , and extended 



the l ocal disc m odel of 
disc. Dwek et al. 



Zheng et al 



(2001) to the whole 



tested a series of models against images 
of the Galactic bulge from the Cosmic Background Explorer 
(COBE) satellite, and found their G2 model to provide one of 
the best fits. This model bar extends from 3-13 kpc and is in- 
clined to the Gala ctic centre line of sight (LOS) at an angle 

Of #bar 



13.4°. 



Gerhard (2002) states that physical mod- 



els can be found for the COBE bar with angles in the range 
15° J; ^bar IS 35°, and many studies assume #b ar ~ 20°. 
However, more recent data from GLIMPSE (Galactic Legacy 
Infrared Mid-Plane Survey Extraordinaire), using the SST, 



support a much larger value of (44 ± 10) ° (Benjam in et al 



2005), while from EROS-2, 



Hamadache et al 



(2006) report 



9bar = (49 ±8)°, which is consistent with original OGLE-I 
results ( Stanek et al ] |l994l) . Hence, predictions are also made 



here for 



Dwek et al 



s E2 model, which has the largest bar an- 



gle of their models: #b ar — 41.3°. 
2.2 Source population 

The expected togle, tmacho and tmoa are to be calcu- 
lated. Therefore for each LOS, the apparent magnitude dis- 
tribution of the model sources must match the observed dis- 
tribution 



Sumi (2004) fitted the /-band stellar distributions 
in 48 OGLE-II Galactic bulge fields with the power-law plus 
Gaussian luminosity function 

(I ~ Wrc) 2 



Mi) =poW pi1 + 



P2 exp 



?rr 2 
Zcr /,RC 



(1) 



where po, pi, pi and o\t,rc are free parameters, and {I)rc 
is measured as described in his paper. The power-law part 
contains red giants and bright main-sequence stars, which 
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m 


6(°) 


Angular separation (° ) 


OGLE 


1.16 


-2.75 




Field 34 


1.35 


-2.40 


0.40 


MACHO 


1.50 


-2.68 




Field 20 


1.68 


-2.47 


0.28 


MOA 


3.0 


-3.8 




Field 36 


3.16 


-3.20 


0.6 



Table 1. Selection of the OGLE-II Galactic bulge fields that are clos- 
est to the lines of sight of the OGLE, MACHO (RCG) and MOA op- 
tical depth measurements. 



lie throughout the bar. The Gaussian component consists of 
RCGs, which in the model here are more concentrated in the 
central part of the bulge, occupying the region 6-10 kpc. (This 
concentration is found to improve the match to the observed 
magnitude distributions, and is not unreasonable, as RCGs are 
older, evolved stars, and hence more likely to exist only in 
mor e densely pop ulated regions). 



Sumi 



1 2004) thus provides, for each of these fields, an 



observed distribution of apparent magnitude. The posi tions of 



Udalski et al 



(2002). 



these fields are listed in table 1 of 1 

The fields closest to the OGLE, MACHO (RCG) and 
MOA lines of sight are selected, as shown in Table [T] The 
MACHO (DIA) LOS is not considered, as explained below. 

As described in §2.3\ the apparent magnitude of each 
source is calculated by first assigning it an absolute mag- 
nitude, and then correcting for its distance. Hence for each 
LOS a separate model distribution of absolute magnitude is 
required that will, with distance corrections, reproduce the ob- 
served distribution of apparent magnitude. Of course in real- 
ity the absolute magnitude distribution should be virtually the 
same for each direction in Table [T] since over these small an- 
gular separations the mass function is expected to vary little. 
Here, the artificial absolute magnitude distributions are only 
used as a means to ensure the model distributions of appar- 
ent magnitude match those observed for each direction. It is 
assumed that the forms of the two distributions are the same, 
i.e. a power-law plus Gaussian. For each of the three OGLE- 
II fields listed in Table Q] an appropriate absolute magnitude 
distribution can easily be generated, by suitably adjusting the 
(extinction-corrected) fitted parameters o f the apparent mag- 



nitude distribution found by 





-^0,min 


-^0,max 


OGLE 


12.1 


15.3 


MACHO 


13.9 


16.2 




(V = 16.37) 


(V = 20.19) 


MOA 


13.6 


20.8 



Table 2. Defined ranges of detectable, extinction-corrected apparent 
magnitudes for OGLE, MACHO and MOA. 



are not given in Sumjj I 12004D . and are provided by T. Sumi, 
priv ate commun i cation ) . 



Sumi et al 



Sumi et al. 



1 2006L table 4) 



(2003, table 5) and[! 
list the extinction-corrected /-band apparent magnitudes of all 
the MOA and OGLE sources used in their respective measure- 
ments of r. The minimum and maximum magnitudes given in 
each case a re taken to define rang es of detectable apparent 



magnitudes. 



Popowski et al 



(2005), in their table 2, provide 
uncorrected l/-band apparent magnitudes. Since these MA- 
CHO sources are all RCGs, their apparent magnitudes are con- 
verted to / - band using the following relation for RCGs from 



Sumi et al 



(2003): 



1= (1.45 ±0.12) (V- I) + 12.7. 



(2) 



MACHO source extinction is then accounted for by simply 
shifting the minimum and maximum MACHO magnitudes by 
the mean A i for the corr esponding OGLE-II field, as given 
in table 3 of ! 



' Sumil d2004l) : for field 20, Ai = 0.951. (Strictly 
speaking the extinction should instead be calculated for each 
source individually, but this is neglected as in the region occu- 
pied by the model RCGs, Ai varies from the mean value by 
< 0.0 5 mag - a negligible amount). Although 



Alcock et al 



Sumi ( 2004). (These parameters 



(2000) also list the apparent magnitudes of the MACHO 
sources used in their DIA measurement, these are in V-band 
and do not consist of only RCGs. Therefore an /-band mag- 
nitude range cannot be reliably defined for the model. The de- 
fined ranges of detectable, extinction-corrected apparent mag- 
nitudes Io for OGLE, MACHO and MOA are given in Table 

El 

2.3 Optical depth 

The expected togle, for example, can now be calculated as 
follows. First a distance is chosen for a given source along 
the OGLE LOS. It is assigned an absolute magnitude using 
the artificial distribution constructed for the nearest OGLE-II 
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«,&)(") 


^obs 

(xlO" 6 ) 


7"model,G2 
(xlO" 6 ) 


r modcl,E2 
(XlO" 6 ) 


OGLE 


(1.16, -2.75) 


9 re+0.57 


2.14 


1.57 


MACHO 


(1.50, -2.68) 


9 i 7 +0.47 
Z - L ' -0.38 


2.19 


1.61 


MOA 


(3.0, -3.8) 


3.36tJi 1 1 (2.59«;r 4 ) 


1.38 


1.01 



Table 3. The expected G2 optical depths agree well with those reported by OGLE and MACHO, but not with MOA's values. (The numbers in 
parentheses are without the DIA correction for disc sources - see text). 



field, #34. The source's apparent magnitude is then calculated 
by accounting for its distance. 

If this apparent magnitude falls within the defined range 
of magnitudes detectable by OGLE, the sour ce is included 
in the calculation of r (using equation (5) of 



Wood & Mao 



20051) . This process is then repeated for many sources. The 



expected tmacho and tmoa are similarly calculated. 

3 RESULTS AND DISCUSSION 
3.1 Model results 

Fig.Q]shows that the model well reproduces the observed dis- 
tributions of apparent magnitude for fields 20, 34 and 36. Ta- 
ble [3] shows that the expected values of togle and tmacho 
also agree well with the observed values, for the G2 model. 
However, the expected tmoa lies ~2.4<j below the reported 
value. As the MOA measurement is sensitive to all sources 
along the LOS, a correction was applied to account for disc 
sources. This is expressed by the /disc term in the r measure- 
ments quoted in <JTJ Such adjustments typically raise r by ~25 
per cent. Note that the model underpredicts tmoa by a much 
greater margin, hence the disagreement cannot be attributed 
to the correction applied by MOA. It therefore appears that 
the discrepancy in the survey measurements cannot be simply 
explained by a dependance on their different flux limits. 

However, there may be other ways in which r depends on 
the source flux. So far the predicted optical depths have been 
calculated by summing over all the source stars whose appar- 
ent magnitudes fall within specified ranges. By repeating this 
process for many small bins of la, t can be predicted as a 
function of Io- This is plotted in Fig. [2] for the OGLE, MA- 
CHO and MOA coordinates (togle los, tmacho los and 
tmoa los, respectively). The detectable magnitude ranges 
given in Table [2] are also shown. For each LOS the absolute 



1 



1 



■ Observed (field 34) 
Model (OGLE LOS) 



I 6 



(mag) 



■ Observed (field 20) 
Model (MACHO LOS) 



14 



I 6 



(mag) 



■ Observed (field 36) 
Model (MOA LOS) 



I 6 



(mag) 



20 



20 



20 



Figure 1. (Top, middle, bottom) panel: apparent magnitude distribu- 
tions, for stars observed in OGLE-II field (34, 20, 36) bv lSumil 120041) . 
and model stars along the (OGLE, MACHO, MOA) LOS. These plots 
are for the G2 model - there is negligible difference with E2. The 
model curves are normalised to the same area as the corresponding 
observed curves. 



5 





3 







' ' 1 ' 


















T 


amplitude, G2 




2 


5 
























G2 model I 


o 


2 















6 














E2 model 




1 


5 










o 
h 












7 t amplitude, E2 




1 







Detectable 














range 











5 









12 14 16 18 20 

l (mag) 



3.0 r 




0.5 L , , , , , , , , , , i , , , l_ 

12 14 16 18 20 

l (mag) 



Figure 2. (Top, middle, bottom) panel: expected (togle LOS> 
TVIACHO LOS. TvlOA LOs) as a function of source apparent mag- 
nitude, for the G2 and E2 models. The detectable magnitude ranges 
given in Tableware shown. In the top panel, the amplitude of the 
T OGLE LOS oscillation (see text) is indicated for both the G2 and E2 
models. 

expected value of r is higher for the G2 bar than the E2, but 
its trend with magnitude is similar. These trends are explained 
as follows. 

r increases rapidly over the range 12 < Jo < 13. Al- 
most all sources of magnitude ~12 will be on the near side 
on the bulge, so as Jo increases fainter and more distant stars, 
with higher optical depths, come into view. For Jo > 15.5, r 
is approximately constant. This is because the power-law part 
of the source magnitude distribution spans a wide range of Jo, 



Hence these stars can be either bright or faint whether they 
are near or far, and thus will show little or no correlation be- 
tween apparent magnitude and distance. So, when calculating 
the average r for a given apparent magnitude, the lower op- 
tical depth of the closer stars is balanced by the higher r of 
those more distant. 

In comparison, the Gaussian (RCG) part of the source 
distribution covers only a very narrow range of absolute mag- 
nitudes. The RCGs' distribution in apparent magnitude will 
be broader, due to variations in their distance, but as they are 
more concentrated in the centre of the bulge, this broaden- 
ing is not great. Therefore the vast majority of RCGs will lie 
within a small range of apparent magnitude, and hence show a 
strong correlation between apparent magnitude and distance. 
At Io ~ 14 we see many RCGs, and they greatly outnumber 
the other sources. Most of the RCGs at this magnitude lie on 
the near side of the bulge, and r is lower. As Jo increases, the 
average distance of the RCGs (and so of all sources) being ob- 
served shifts towards the far side of the bulge, and r increases. 
As Jo becomes fainter still, > 15, we see fewer and fewer 
RCGs, and the average distance of all the observed sources 
moves back towards the centre of the bulge, where it then re- 
mains, and t becomes approximately constant. The amplitude 
of this oscillation in r (hereafter the r amplitude) caused by 
the RCGs along the OGLE LOS is indicated in the top panel 

ofFig.m 

This strong correlation displayed by the RCGs is illus- 
trated as follows. In Figs. [3] and [4] contours are plotted of 
source counts as a function of distance and apparent magni- 
tude along each LOS, for the G2 and E2 models, respectively. 
The two components of the source population are clearly dis- 
tinguishable. For the G2 model, most of the RCGs appear as a 
narrow diagonal line in the range 14.0 < Jo < 15.0. For E2, 
this region is wider (for a given distance) and shallower. These 
differences in shape are primarily due to the different bar an- 
gles of the G2 and E2 models (13.4° and 41.3°, respectively). 
The red giants and other stars form a smoother background, 
with a steep increase in numbers, and a broadening in distance, 
as Jo increases. In the top panels of Figs.[3]and[4] the vertical 
dotted lines indicate slices of this distribution - for the OGLE 
LOS - that are shown in Figs.[5]and[6] respectively. Finally, the 
top panel of Fig.|7]gives the average distance of model OGLE 
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Figure 3. (Top, middle, bottom) panel: G2 model source counts as 
a function of distance and apparent magnitude, for the (OGLE, MA- 
CHO, MOA) coordinates. The RCG component is clearly visible (see 
text). The vertical dotted lines (top panel) correspond to the slices 
shown in Fig.[5](see text). The normalisation is arbitrary. Contour lev- 
els are at 0.5, 1.0, 2.0, 5.0, 10.0, 20.0, 50.0 and 100.0. 



LOS sources as a function of 7q- The slice magnitudes are in- 
dicated, and the bottom panel shows how they intersect the 
togle los trend from Fig. [2] Note that the trends of average 
source distance and optical depth with Jo are almost identical, 
as would be expected for the reasons given above. 

The expected oscillation in r caused by the RCGs is 
clearly significant. For example, Fig. [2] shows that the t am- 
plitude along the OGLE LOS is ~1.2 x 1CT 6 . This is a devia- 
tion of ~±30 per cent from the approximately constant optical 
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Figure 4. Model source counts as a function of distance and apparent 
magnitude. Same as Fig. [5] but for the E2 model, with the vertical 
dotted lines corresponding to the slices shown in Fig.[6](see text). The 
normalisation is arbitrary, but consistent with Fig. [3] 



depth at fainter magnitudes (hereafter rfl at ), where far fewer 
RCGS are seen: Tfj at w 2.1 X 10 -6 . For comparison, OGLE's 
measured value of 2.55_to4g x 10 ~ 6 has an uncertainty of 
only ~±20 per cent, so the predicted oscillation ought to be 
detectable if enough sources are observed at the correct mag- 
nitudes. 
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Figure 5. Model source counts as a function of distance (OGLE LOS, G2 model), for selected apparent magnitudes Io (with Iq bin widths of 0.1 
mag). These magnitudes are indicated in the panels, and correspond to the slices indicated in Fig.|5](top panel). 
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Figure 6. Model source counts as a function of distance. Same as Fig. [5] but for the E2 model, with the selected magnitudes corresponding to the 
slices indicated in Fig.[4](top panel). 
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Figure 7. Top panel: average distance of model OGLE LOS sources 
as a function of magnitude. Bottom panel: togle LOS as a function 
of magnitude (same as top panel of Fig. [2). The vertical dotted lines 
correspond to the slice magnitudes indicated in Figs.|3j(6] 





(/,£>) (°) 


#bar,G2 (°) 


#bar,E2 (°) 


OGLE 


(1.16, -2.75) 


14.8 


19.1 


MACHO 


(1.50, -2.68) 


24.6 


31.0 



Table 4. la upper limits on #bar. from combining the expected op- 
tical depths with those measured by OGLE and MACHO (see text). 
Values are shown for the G2 and E2 bar models. 



If the given bar angles of the G2 and E2 models are now 
varied, model-independent trends with f)bar may be revealed. 
Fig. [8] shows that the expected amplitudes of togle los, 
tmacho los and tmoa los all decrease with increasing bar 
angle. This provides a potential constraint on #bar, should the 
expected r amplitude be observed and its magnitude accu- 
rately measured. Another constraint is shown by Fig. [9] where 
the expected optical depths togle, tmacho and tmoa dis- 
play a similar dependance on (9bar- The corresponding ob- 
served values are overplotted, with their la uncertainties, and 
from the intersections with the predicted OGLE and MACHO 
curves, la upper limits on #bar are obtained. (There is no in- 
tersection between the predicted and observed tmoa). These 
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Figure 8. (Top, middle, bottom) panel: oscillation amplitude of ex- 
pected (togle LOS. TMACHO LOS, tmoa LOS) as a function of 
"bar, f° r me and E2 models. 



limits are given in Table|4] Note that they exclude the large bar 
angle of the E2 model (#b ar — 41.3°), as well as those from 
GLIMPSE (#bar = (44±10)°) and EROS ((9 b ar = (49±8)°). 

3.2 Comparison with EROS data 



The EROS-2 survey dHamadache et al 



2006) has found the 



largest sample of clump-giant events so far, 120, compared 
with 32 and 62 for the latest OGLE and MACHO surveys, re- 
spectively. This sample may be sufficient to enable a useful 
comparison of the predicted optical depth trends with obser- 



10 Wood 




3.0 [ 



40 60 
Bar angle (°) 




20 



40 60 
Bar angle (°) 




40 60 
Bar angle (°) 



30 



Figure 9. (Top, middle, bottom) panel: expected (togle, tmacho. 
tmoa) as a function of 6>bar. f° r me G2 and E2 models. The horizon- 
tal lines show the value measured by (OGLE, MACHO, MOA) with 
its ltT uncertainties. 
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Figure 10. (Top, middle, bottom) panel: expected (togle LOS. 
TMACHO LOS. TMOA LOs) as a function of source apparent mag- 
nitude, with and without stellar flux smearing of 20 per cent rms, as 
indicated. The curves with no smearing are the same as in Fig. [2] 



vational data. Since the EROS-2 survey (like all microlensing 
surveys) was conducted across many fields rather than for a 
specific LOS, any observed oscillation in r similar to the pre- 
diction would be somewhat smoothed out. However, this ef- 
fect should not be strong, as the predicted trend is similar for 
different lines of sight towards the bulge (as shown in Fig.O, 
and in any case such effects are accounted for as described 
below. 

To make the comparison, unpublished EROS-2 data have 
been supplied by J. Rich et al. (private communication). They 



find, from studies with artificial stars, that clump giant fluxes 
are smeared by ~20 per cent rms, which does not affect their 
optical depth calculations averaged over the whole clump, but 
will of course reduce the r amplitude by smoothing out the 
predicted oscillation. This flux smearing effect is added to all 
of the model stars, using a Gaussian with a = 0.2. It is found 
that although the r amplitude is indeed reduced, the oscillation 
is still clear, as shown in Fig. 1101 

EROS stars were observed in two non-standard bands, 
-Reros and Beros, where .Reros = /ogle- Hence only the 
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J?eros data are considered, to avoid any c omplications aris 



Hamadache et al 



ing fro m using different magnitude scales. 
(2006) divided each of their 66 bulge fields into 32 subfields. 
For each subfield, they modelled the stellar density in colour- 
magnitude space as a power-law plus Gaussian, in order to 
find the magnitude, J?eros, clump, and colour of the clump 
centre. Hence smoothing effects (as mentioned above) and ex- 
tinction can now both be easily accounted for here, by taking 
the clump centres as reference points: instead of simply find- 
ing r as a function of J?eros, it is found as a function of mag- 
nitude relative to the clump centres, J?eros — ^eros, clump- 
This method 'lines up' all the clump centres in each subfield. 
For comparison, the predicted r is similarly found as a func- 
tion of J - Jo, dump. For the (OGLE, MACHO, MOA) LOS, 
the Gaussian component shown in Fig.Qjpeaks at Jo, clump = 
(14.75, 14.65, 14.55) ± 0.05 (for bin widths of 0.1 mag). For 
the E2 model, the magnitudes are slightly brighter: Jo, c ium P 
= (14.65, 14.55, 14.45) ± 0.05. 

teros is now foun d as a function of source magnitude 
using equation (13) from 



Hamadache et al 



(2006) 



f Si=l te,i/e(tB,i 



2uo (max) 



(3) 



where each event i has a time-scale tE,% and detection effi- 
ciency e(tE,i), each monitored star j is observed for a time Tj, 
the total numbers of events and stars are iVcv and N* , respec- 
tively, and the maximum impact parameter no (max) = 0.75. 
This calculation is implemented with a separate summation 
for each magnitude bin (and taking full account of the differ- 
ent detection efficiencies for each EROS-2 field). The uncer - 
tainty is also determined following 



Hamadache et al 



(2006), 



who added in quadrature a 5 per cent systematic part - due 
to blendin g effects - and a large r statistical part, estimated ac- 



Han& Gould (1995) 



cording to 

Fig. QTJ shows teros as a function of J?eros — 
J^eros, clump- The observed trend is now compared with the 
flux-smeared predicted trends, using a \ 2 test, to see if the for- 
mer is better fitted by an oscillating or constant optical depth. 
Since it is fitting an oscillating form that is of interest, the 
width of the oscillation and the absolute values of the r ampli- 
tude and Tfl a t are allowed to vary as free parameters in the fit: 
the model trend may be stretched along the magnitude axis, 
and both shifted and stretched along the r axis - though not 



stretched in the former case by more than an arbitrary limit of 
50 per cent, since Fig.[2]shows that the width of the oscillation 
does not vary much with direction. Finally, to allow for other 
slight changes in the shape of the oscillation with direction 
and bulge model, all six of the predicted trends (three lines of 
sight, G2 and E2 models) are tested. These x 2 values are com- 
pared with that for a freely-fitted constant optical depth. The 
results are shown in Table [5] Also indicated are the chance 
probabilities p that \ 2 would be greater than or equal to the 
given values. 

An oscillating r appears to provide a better fit to the 
data than a constant optical depth. There is a mostly negli- 
gible change in \ 2 with direction, and a small but insignifi- 
cant preference for the E2 model. Fig. QTJ shows the best-fit 
oscillation, with the MACHO E2 trend: r na t has been shifted 
by -0.05 x 10~ 6 , and the curve has been stretched by fac- 
tors of 1.50 and 1.60 along the magnitude and r axes, respec- 
tively. (Note that the factor of 1.50 is at the (arbitrary) 50 per 
cent limit. The fits improve with further magnitude stretching, 
the best possible fit being for a factor of 2.50 (MACHO G2 
model), with \ 2 = 1-44, but this is well beyond the limit and 
is ignored). However, the significance of the % 2 preference for 
an oscillating r, rather than a constant value, is low. A reason- 
able magnitude binning gives only eight data points. Whereas 
the oscillating r fit has three free parameters and five degrees 
of freedom, the constant r fit has of course just one free pa- 
rameter, and seven degrees of freedom. The constant r fit is 
not significantly discrepant to the data. 

It is possible that the EROS-2 event detection efficiency 
may be a function of magnitude. However, it would not be 
a strong dependance, and any variation would be smooth. It 
could not therefore hide any real oscillation of r, or generate 
a false one (J. Rich, private communication). It appears that 
the available data are still insufficient to accurately determine 
the dependance of the optical depth on source apparent mag- 
nitude. 

However, there is another, simpler way to look for signs 
of the predicted r oscillation in the survey data. If r is indeed 
higher on the faint side of the clump centres, then more of the 
observed events should also be on the faint side: a plot of event 
counts, as a function of J?eros — J^eros, clump, should be no- 
ticeably skewed towards the faint side. Fig. Q2] shows that this 
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Figure 11. Observed teros> and fitted model tmacho LOS. as 
functions of magnitude relative to the clump centre, Reros — 
^EROS, clump and I - 1 0> c i ump , respectively (see text). The latter 
magnitude scale is stretched relative to the former by the fitting factor. 
Also shown is the fitted rg at . 
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Figure 12. EROS event counts as a function of magnitude relative to 
the clump centre, Reros — Reros, clump- The error bars represent 
Poisson uncertainties. The distribution is skewed towards the faint side 
of the magnitude axis; the trend is fitted with a Gaussian in order to 
test the significance of the skew (see text). The vertical dotted line 
indicates Reros — Reros, clump = 0, and the vertical dashed line 
indicates the centre of the Gaussian. 





G2 


E2 


OGLE 


2.07 


1.84 


MACHO 


1.91 


1.83 


MOA 


1.99 


1.93 


Average 


1.99 


1.87 




(p > 0.75) 


(p > 0.75) 


Constant r 


4.24 




(p > 0.75) 


Ax 2 


2.25 


2.37 



Table 5. \ 2 values from fitting the observed teros as a function of 
Reros ~ fidump with the predicted (oscillating) trends, for differ- 
ent lines of sight and bulge models as indicated, and with a constant 
optical depth. Also indicated are the chance probabilities p that \ 2 
would be greater than or equal to the given values. An oscillating r 
provides a better fit, but at a low significance (see text). 



is so. The ratio of events with -Reros — R^eros, clump fainter 
than zero to those brighter than zero (hereafter the skewness 
ratio) is 1.14. The significance of this skew is tested as fol- 
lows. The plot of event counts is fitted with a Gaussian, as 
shown. This is of course an imperfect fit, but gives an ap- 
proximate measure of the dispersion. Then 120 points (for the 
120 EROS events) are randomly distributed on the magnitude 
axis, according to a Gaussian with the fitted dispersion, and 
the skewness ratio is calculated. For one million such Monte 
Carlo cases, only 21 per cent have a skewness ratio greater 
than the observed value. Thus the observed skew is sugges- 
tive, but not highly significant. 



4 SUMMARY AND CONCLUSIONS 

It does not appear that the discrepancy in optical depth mea- 
surements between the RCG and all-star analyses can be ex- 
plained by a dependence of the lensing surveys on their flux 
limits. The model reproduces the OGLE and MACHO val- 
ues based on RCGs, but underpredicts MOA's all-star value 
by <-~2.4o\ Another potentia l explanation for the discrepancy 
is blending. 



Sumi et al 



(2006) found ~38 per cent of OGLE- 
II events with apparent RCG sources were really due to faint 
stars blended with a bright companion. However, they also 
showed that blending has little effect on estimates of r due 
to parti al cancellation of its d iffe rent effects, a point als o 
made by 



Popow ski et al. 



Sumi et al 



Hamadache et al 



(2006) 



( 20051) and|] 

(2006) also state that the DIA method is less sensi- 
tive to the systematics of blending in crowded fields. Though 
it is of course possible that MOA's value may yet be low- 
ered with more data, it is supported by MACHO's earlier DIA 
value. 

r is expected to be generally constant as a function of 
source apparent magnitude for Iq > 13.0, except in the range 
13.5 < Io < 15.5, where many RCGs are detected. These 
stars dominate the source counts at such magnitudes, and show 
a strong correlation between distance and apparent magni- 
tude, causing a significant oscillation in r. The amplitude of 
this oscillation is found to decrease with increasing bar angle, 
providing a potential constraint on #bar. A further constraint 
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comes from a similar dependance of r with 8^^: combining 
the predicted trends with the measured values provides la 
upper limits, which exclude the large bar angles reported by 
GLIMPSE and EROS. 

From the EROS-2 survey, teros has been found as a 
function of source apparent magnitude. The predicted oscilla- 
tion is not only consistent with the observed trend, but pro- 
vides a better \ z fit than a constant optical depth, though 
the significance of this preference is low due to insufficient 
data. However, a further sign comes from EROS event counts, 
which show a clear skew towards fainter magnitudes. With 
ongoing surveys detecting increasing numbers of RCG events 
(and ~500 yr _1 of all kinds), it should soon be possible to 
make a more useful and definite comparison. 
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